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ABSTRACT 

o ■ 

^ , The X-ray source IE 1207.4—5209 is a compact central object in the 

pH . G296. 5+10.0 supernova remnant. Its spin period of 424 ms, discovered with 

^ . the Chandra X-ray Observatory, suggests that it is a neutron star. The X-ray 

^ , spectrum of this radio-quiet pulsar shows at least two absorption lines, first spec- 

^ ! tral features discovered in radiation from an isolated neutron star. Here we report 

^ ! the results of timing analysis of Chandra and XMM-Newton observations of this 

^ ! source showing a non-monotonous behavior of its period. We discuss three hy- 

0\ ', potheses which may explain the observational result. The first one assumes that 

^ ! IE 1207.4—5209 is a glitching pulsar, with frequency jumps of A/ > 5 /iHz oc- 

^ ! curring every 1-2 years. The second hypothesis explains the deviations from a 

Q ! steady spin-down as due to accretion, with accretion rate varying from ~ 10^^ to 

r-| ! > 10^^ g s~^, from a disk possibly formed from ejecta produced in the supernova 

^! explosion. Finally, the period variations could be explained assuming that the 

O I pulsar is in a wide binary system with a long period, Porb ~ 0.2-6 yr, and a 

^ ', low-mass companion, M2 < O.SMq. 

• 1— I ! Subject headings: pulsars: individual (IE 1207.4—5209) — stars: neutron — 

^ ! supernovae: individual (PKS 1209—51/52) — X-rays: stars 



1. Introduction 

X-ray observations of supernova remnants (SNRs) have revealed a number of radio- 
quiet central compact objects (CCOs — see Pavlov et al. 2002a, 2004 for a review) whose 
nature remains enigmatic. They are characterized by soft, apparently thermal. X-ray spectra 
and a lack of manifestations of pulsar activity (e.g., radio and/or 7-ray emission, compact 
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synchrotron nebulae). Most likely, CCOs are neutron stars (NSs) formed in supernova ex- 
plosions. 

One of the best investigated CCOs is IE 1207.4-5209 (1E1207 hereafter) in the G296.5+10.0 
(= PKS 1209-51/52) SNR. It was discovered by Helfand & Becker (1984) with the Einstein 
observatory. Mereghetti, Bignami & Caraveo (1996) and Vasisht et al. (1997) interpreted the 
ROSAT and ASCA spectra of 1E1207 as blackbody emission of a (redshifted) temperature 
kT'^ ~ 0.25 keV from an area with radius R°° ^1.5 km ai d = 2 kpc. Mereghetti et al. 
(1996) put upper limits of ~ 0.1 mJy for 4.8 GHz radio flux, 10^^ photons cm^^ s^^ for 
E > 100 MeV 7-rays, and V > 25 for an optical counterpart. Zavlin, Pavlov & Triimper 
(1998) showed that the ROSAT and ASCA spectra are consistent with hydrogen or helium 
NS atmosphere models. For a NS of mass 1.4 and radius 10 km, they obtained a NS 
surface temperature kT^ = kT {1 + z)~^ — 0.12-0.16 keV, where z is the gravitational 
redshift at the NS surface, and a distance 1.6-3.3 kpc, compatible with the distance to the 
SNR, d = 2Atli kpc (Giacani et al. 2000). 

First Chandra observation of 1E1207 in January 2000 with a 30 ks exposure allowed 
us to detect a period Pi — 0.4241296 s (frequency /i = 2.357770 Hz), which proved that 
the source is indeed a NS, with a period typical for an ordinary radio pulsar (Zavlin et al. 

2000; Paper I hereafter). The pulsar was again observed with Chandra for 30 ks in January 
2002. This observation showed a longer period, P2 = 0.4241309 s (/s = 2.357763 Hz), 
corresponding to the period derivative P ^ 2 x lO^^"^ s s"^ (Pavlov et al. 2002b; Paper II 
hereafter). Such a period derivative implies a characteristic age, = P/{2P) ~ 300 kyr, 
much larger than the 3-20 kyr age of the SNR (Roger et al. 1988), which suggests that the 
pulsar was born with a period only slightly faster than its current value. The conventional 
magnetic field, S ~ 3 x 10^^ G, and the spin-down energy loss rate, £^ ~ 1 x 10^^ erg s~^, 
inferred from the P, P values, are typical for radio pulsars. 

Even more surprising finding from the two Chandra observations was the discovery of 
absorption fines, first lines detected in the spectrum of an isolated NS (Sanwal et al. 2002a). 
Two firmly detected lines, at 0.7 keV and 1.4 keV, could be interpreted as absorption fines 
of once-ionized helium in a very strong magnetic field, about 1.5 x 10^^ G, which requires a 
gravitational redshift z = 1.12-1.23 (Sanwal et al. 2002a). Another interpretation, that the 
lines could be associated with transitions in He-like oxygen ions in a magnetic field B ~ 10^^ 
G, was discussed by Hailey & Mori (2002). A possible third line, at about 2 keV (Sanwal 
et al. 2002a), unfortunately coincides with the Ir M line from the telescope mirror coating, 
where the calibration is inaccurate. If this line is due to the source, then the three lines might 
be the fundamental and two harmonics of the electron cyclotron absorption in a magnetic 
field ~ 0.6 (l+z) X 10^^ G. Sanwal et al. (2002a) discussed this possibility and concluded that 
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it is difficult to reconcile this interpretation with the much higher magnetic field estimated 
from P and P, and with the very low oscillator strengths of cyclotron harmonics. 

First observation of 1E1207 with XMM-Newton (27 ks exposure) was performed in 
December 2001. Mereghetti et al. (2002) confirmed the absorption lines and reported on 
their pulse phase variations. The period value was in agreement with that measured in the 
second Chandra observation, taken 13 days later. 

1E1207 was again observed with XMM-Newton in August 2002 (258 ks of total expo- 
sure). Based on this very deep observation, Bignami et al. (2003) reported positive detection 
of the 2.1 keV absorption line and marginal detection of a 2.8 keV line in the source spectrum. 
The detection of three (possibly four) evenly spaced lines argues for the above-mentioned 
cyclotron interpretation, although it remains unclear how the cyclotron harmonics can form 
in the relatively cold plasma with the low magnetic field. Timing analysis of these data (de- 
scribed below) provides a highly accurate period that is shorter than the period measured in 
the first XMM-Newton observation, indicating that the pulsar is not spinning down steadily. 

Finally, Chandra observed this puzzling pulsar two times in June 2003 for 280 ks of total 
exposure. These observations were primarily designed to measure the phase-dependent spec- 
trum with a high energy resolution. However, they also provided timing information that we 
use to assess the evolution of the pulsar's period. We describe the XMM-Newton and Chan- 
dra observations and present results of our timing analysis in § 2. Possible interpretations 
are discussed in § 3. The results are summarized in § 4. 

2. Observations and timing analysis 

2.1. XMM-Newton observation of December 2001 

The timing analysis of the first (2000 January 6-7) and second (2002 January 5-6) 
Chandra observations has been described in detail in Papers I and II. First timing results 
for the XMM-Newton observation of December 2001 have been presented by Mereghetti et 
al. (2002). To evaluate the most plausible frequency and its uncertainty in a statistically 
rigorous way, and for the sake of uniformity, we reanalyzed these data with the same approach 
as in Papers I and II, using the method of Gregory & Loredo (1996) based on the Bayesian 
formalism. While this method yields frequency estimates consistent with those given by other 
techniques (e. g., the simple epoch-folding and Rayleigh tests), it has several important 
advantages. It is free of any assumption on pulse shape, that results in more accurate 
determination of signal frequency and its uncertainty. The Bayesian approach implements 
the phase-averaged epoch-folding algorithm to calculate the frequency-dependent odds-ratio 
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0{f) that allows one to find the corresponding probability distribution p{f) oc 0{f), which 
in turn is very useful for the interpretation of the results (see Gregory & Loredo 1996 and 
Papers I and II for more details). In the following, we count all frequencies from a reference 
value /rcf = 2.3577 Hz. For each of the epochs, we measure the mean frequency /mean with 
the standard deviation a and provide them in the form / = /mean ± o" (e.g., ff^ = 69.9 ± 1.3 
/iHz and /|'^ = 62.5±3.7 /^Hz for the first and second Chandra observations — see Paper II). 
We also measure the median frequency /med with the 68% and 90% uncertainties, cr^ and 
(7^, below and above /med) and provide them in the form / = /med(~cr68' +"^68' ~'^90' +"^90) 
[e.g., /f = 69.9 (-1.4, +1.2; -2.5, +1.8) //Hz, f^^ = 62.5 (-3.2, +4.5; -4.2, +8.0) /iRz]. 

First observation of 1E1207 with XMM-Newton was performed on 2001 December 23- 
24 (orbit 374). For the timing analysis, we used the data taken with the European Photon 
Imaging Camera based on the 'p-n' CCD technology (EPIC-pn) operated in Small Win- 
dow mode, which provides a 4'4 x 4'4 sky image and a 5.7 ms time resolution. The total 
observation time span was 27.0 ks, corresponding to an effective exposure of 18.9 ks (be- 
cause of about 30% time loss during the CCD readout). The data were processed with the 
most recent version of the "oal" library (v. 3.108)^. We used 26 778 counts extracted from 
a 30"-radius circle centered at the source position in the 0.2-5.0 keV energy range. We 
estimated that about 88% of those counts belong to the source. The photon arrival times 
were transformed to the Solar System Barycenter with the "barycen" task. The test (see 
Papers I and II for details) gives the most significant peak, Zf — 90.9, at the frequency 
of 63.1 //Hz, very close to 63.0 ± 2.6 /(Hz found by Mereghetti et al. (2002) from the 
test assuming that the signal has a sinusoidal shape. Implementing the Bayesian approach 
results in the probability density distribution p{f) shown in Figure 1. The mean and median 
frequencies of the corresponding probability distribution are /™™ = 59.2 ±1.3 /xHz and 
jxmm _ 2 (—1.4^ +0.9; —2.5, +1.5) /iHz. These values are in agreement with the frequen- 
cies obtained from the Zf and tests. They do not show statistically significant differences 
with the frequency (see above) found from the second Chandra observation 13 days later. 



2.2. XMM-Newton observations of August 2002 

XMM-Newton observed 1E1207 again on 2002 August 4-5 (orbit 486) and 6-7 (orbit 
487) with the same instrumental setup as in the observation on orbit 374. We excluded 
intervals of strong background ffares at the beginning and at the end of each observation 
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and used uninterrupted spans of 107.0 and 111.0 ks from the data taken on orbits 486 
and 487, respectively (the total effective exposure of EPIC-pn is 152.6 ks). The time gap 
between these two intervals is 73.0 ks. Applying the same extraction radius and energy 
range as for the orbit 374 observation, we obtained 108 945 and 113 720 counts from the two 
data sets for the timing analysis (background contamination is about 11%). The frequency 
distributions of the odds-ratios give mean frequencies /2™™ = 63.31 ± 0.48 and = 
64.18 ± 0.33 /iHz, and median frequencies f^^"^ = 63.19 (-0.26, +0.78; -0.58, +0.14) /xHz 
and f^^"^ = 64.30 (-0.53, +1.08; -0.71, +0.39) //Hz for the orbits 486 and 487, respectively. 
These estimates yield the frequency shift between the two orbits, A = +0.88 ± 0.59 //Hz 
(calculated from the probability distribution for the frequency difference — see Paper II). 
The probability that the shift is positive is P(A > 0) = 0.932. 

This shift, albeit of low statistical significance, indicates a positive frequency derivative, 
/ ~ 10~^^-10^^^ s~^, in the joint data set (Tgpan = 291.0 ks). To examine the possible 
frequency change, we calculate the odds-ratio 0(/, /) and the Raylcigh statistic Zf{f, f) on 
a two-dimensional (/, /) grid, assuming a linear change of the spin frequency within this 
observation and choosing the middle of the total observation time span as the reference epoch 
to minimize the /-/ correlation. The Rayleigh test gives Zf^^^^ = 448.1 at / = 63.41 /iHz, 
/ — +1.7 X 10^^^ s^^. To estimate the mean and median parameters and their uncertainties, 
we used the probability density distribution g(/, /) = A 0(/, /), where A is the normahzation 

_ —.xmm 

constant, and obtained /f™"^ = 63.49 ± 0.08 //Hz, = (+2.1 ± 2.0) x 10"^^ 3-2^ 

= 63.51 (-0.09, +0.06; -0.16, +0.10) /iHz , 
yxmm ^ +2.0 (-1.8, +2.2; -2.4, +2.6) X 10"^^ s-2. (1) 

These estimates suggest that the frequency was indeed increasing during that observation 
j-p^yxmm ^ g-j _ 0.908], but the Uncertainties are too high to conclude this firmly. The 
probability density distribution for frequency, p{f) — J q{f, f) df, is shown in Figure 1. 

The probability distribution for the frequency difference between the XMM-Newton 
observations of August 2002 and December 2001 is shown in the upper panel of Figure 2. The 
mean and median values for the difference are A = +4.30 + 1.28 /tHz and A = +4.26 (—0.81, 
+1.08; —1.69, +2.31) /(Hz. The probability that the pulsar has spun up in about 7.5 months 
since the first XMM-Newton observation is P(A > 0) = 0.992. If the frequency were 
increasing monotonously during this period (1.9 x 10^ s), then its derivative would be in the 
range / = (1.4 - 3.4) x 10~^^ s~^ (at a 90% confidence level). 
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2.3. Chandra observations of June 2003 

1E1207 was observed with the Low Energy Transmission Grating in combination with 
the Advanced CCD Imaging Spectrometer operated in Continuous Clocking mode on 2003 
June 10-12 (155.7 ks) and June 18-19 (115.2 ks). This observational mode provides a 2.9 ms 
time resolution by means of sacrificing spatial resolution in one dimension. For the timing 
analysis, we used 11909 and 8804 counts (for the first and second data sets, respectively) 
extracted from segments of a 7-pixel (3'/44) width in the zero-order images, in the energy 
range of 0.2-5.0 keV. The times of arrival were corrected for the dither and the Science 
Instrument Module motion as described in Paper I and transformed to the Solar System 
Barycenter using the "axBary" tool of the CIAO package^. 

The mean and median frequencies, as given by odds- ratios, are /g^ = 64.11 ± 0.96, 
fl^ = 62.47 ± 1.36 /iHz, and = 64.35 (-1.40, +0.64; -1.93, +1.00) /iHz, f^^^ = 62.09 
(-0.59, +1.78; -1.04, +3.23) fiHz (the subscripts '3a' and '3b' are related to the observations 
performed on June 10-12 and 18-19, respectively). The frequencies of the Zf peaks, 64.0 
/iHz (^f max — 27.5) and 63.5 /iHz (2'^ max — 25.4), are within the uncertainties of the mean 
and median frequencies. The mean and median values of the frequency difference between 
the two observations are A = -1.67 ± 1.62 //Hz and A = -1.98 (-1.10, +1.84) //Hz. 
The probability that the pulsar has spun down in a week between the two observations is 
P(A < 0) = 0.853. 

We also performed the timing analysis of the combined data set (Tgpan = 799.9 ks) on 
a two-dimensional grid of / and /, using the same approach as in §2.2. The maximum 
value of the Rayleigh statistic is Zl^^ = 52.7, at / = 63.61 //Hz, / = -1.0 x lO^^^ 
s~^. Because of the large time gap, Tgap = 529.1 ks, the phase coherence between the two 

data sets was lost, that resulted in a number of peaks (frequency aliases) in the frequency 
dependences of O and Zf, separated by ~ 1.4-1.5 /iHz. Three most significant peaks (we 
denote them as A, B, and C) are seen in the probability density distribution p{f) in the 
bottom panel of Figure 1. The mean and median parameters as estimated from the the multi- 

-.ch 

peak two-dimensional distribution are ff^ = 63.37+0.74 /xHz, /g = (—3.0 + 2.1) x 10^^^ s^^, 
and /l^ = 63.58 (-0.10, +0.07; -1.37, +0.42) //Hz, /|'^ = -3.5 (-1.6, +2.5; -2.7, +3.9) x 
10~-'^^s~^. Considering (/, /) domains around the peaks A, B and C separately results in 
the following mean and median frequencies: /a — 62.10 + 0.05 /tHz, fs — 63.59 + 0.05 /(Hz, 
fc = 64.96 ± 0.07 /(Hz, 

fA = 62.09 (-0.03, +0.04; -0.08, +0.10) //Hz, 
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/b = 63.60 (-0.06, +0.04; -0.16, +0.11) //Hz, 

fc = 65.02 (-0.06, +0.09; -0.12, +0.17) //Hz, (2) 



and frequency derivatives: /a-u — 


-1.6 ± 


1-4, /b,-12 


= -3.6 ± 1.^ 


/c,-i2 = +0.4±1.6. 


/a,-12 = 


-1.4(- 


-1.2, +0.9; - 


-2.0, +1.7), 




/b,-12 = 


-3.8 (- 


-1.2, +1.6; - 


-1.8, +2.6), 




fc -12 — 


+0.5 (- 


-1.5, +1.1; - 


-2.5, +1.7), 


(3) 



where /_i2 = //(lO^^^ s~^). The relative contributions of the peaks A, B, C into the prob- 
abihty distribution (i.e., the probabihties that the true frequency and its derivative are 
associated with a given peak) arc P4 = 0.207, Pb = 0.734, Pc = 0.059. Therefore, ahhough 
the parameters related to peak B are more probable, none of the two other parameter sets 
can be ruled out on statistical grounds. 

The probability distribution of the frequency difference between June 2003 and August 
2002, plotted in the lower panel of Figure 2, also has 3 discernible peaks. The mean frequency 
shift in the 10 month period is only A = +0.06 ± 0.61 /iHz; the median frequency difference 
is A = +0.09 (-0.12,+0.10; -1.46, +0.14) //Hz. If we assume that peak A (or B or C) 
corresponds to the true frequency in June 2003, then the shifts are A^ = —1.36 + 0.10 
/(Hz and A^ = -1.37 (-0.07,+0.09; -0.10,+0.16) /tHz, or A^ = +0.10 ± 0.07 /(Hz and 
Ab = +0.11 (-0.05,+0.07; -0.09,+0.12) /(Hz, or Ac = +1.47 + 0.09 /(Hz and Ac = +1.48 
(-0.08, +0.07; -0.14, +0.18) /(Hz. Assuming a linear time dependence of frequency during 
this period (2.7 x 10'' s), we can constrain the frequency derivative, —5.5 x 10~^^ < / < 
+6.0 X 10"^^ s-2, at a 90% confidence level. 



3. Discussion 

It is easy to check that the time dependence of the pulsation frequency in the 3.45 yr 
interval cannot be satisfactorily fitted by a straight line, f{t) = fo + fit — to)) for any of 
the three possible frequency values obtained in the third Chandra observation. The best fit, 
with /o = 68.0 /xHz, / = -5.2 x 10"^'' s^^ (to = 51500.0 MJD), is obtained assuming the 
correct frequency in June 2003 is given by peak A (see the dotted line in Fig. 3); the fit 
corresponds to xLin — 22.6 (for 3 degrees of freedom) and can be rejected at a 4.1(7 level (if 
one assumes the frequency of peak B for the third Chandra observation, then the rejection 
level is 5.6a). This means that 1E1207 is not spinning down steadily, as most radio pulsars 
do. We discuss possible explanations below. 
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3.1. A glitching pulsar? 

Some radio pulsars and AXPs occasionally show sudden increases in pulsation frequency, 
commonly known as "glitches", with various patterns of post-glitch behavior (e.g., Lyne, 
Shemar & Graham-Smith 2000; Gavriil & Kaspi 2002). The relative frequency jumps. A///, 
vary from ~ 10~^ to ~ 10~^ in different pulsars. Obviously, the non-monotonous behavior 
of f{t), inferred from the above-described timing of 1E1207, may suggest that the pulsar 
experienced a number of glitches during the 3.45 yr time span. Assuming the glitches are 
not associated with a substantial change of frequency derivative /, at least two glitches are 
required, between January and August 2002 (212 days interval), and between August 2002 
and June 2003 (313 days interval). The values of the frequency jumps depend on the value of 
/ assumed, which in turn depends on the cumulative frequency change in glitches that might 
occur in the 716 days interval between January 2000 and December 2001. For instance, if 
there were no glitches in that interval, then / ~ —1.6 x 10~^^ s~^, and the two frequency 
jumps, of about 7 //Hz and 5 //Hz, are required to fit the data (the latter value assumes that 
the correct frequency in 2003 June 10-19 is given by the central peak B of the probability 
distribution — see Fig. 3). The larger of these jumps (A/// ~ 3 x 10~^), is similar to the 
strong glitches observed in the Vela pulsar (Dodson, McGuUoch & Lewis 2002) and the AXP 
IE 2259+586 (Kaspi et al. 2003), and it is a factor of 5 smaller than the giant glitch recently 
observed in PSR J1806-2125 (Hobbs et al. 2002). 

Given that at least two glitches are required in the 525 day interval, it seems reasonable 
to assume that the pulsar was also glitching in the 717 day interval between the first two 
observations, which would correspond to a larger |/| and stronger glitches. For instance, if the 
cumulative frequency increase due to glitches was 30 fiHz in 717 days, then / ^ —6.6 x 10~^^ 
s^'^, and the two frequency jumps (in 525 days after January 2002) are about 16 //Hz and 18 
/(Hz (A/// ~ 7 X 10~^). Although the relative frequency jumps of the separate glitches are 
smaller than for the largest ghtch observed (Hobbs et al. 2002), the integrated amphtude 
of the glitches is a factor of 30-60 greater then the typical value, J2^f ~ 0.02|/|T {T 
is the total time span, equal to 3.45 yr in our case) observed for glitching radio pulsars 
(Lyne et al. 2000). Therefore, although we cannot formally rule out the possibility that 
1E1207 experiences glitches of A/ > 5 /iHz, with a characteristic rate of ~0.5-l glitches 
per year, such a hypothesis implies that the nature of the glitches in this unusual object is 
different from that observed in radio pulsars. In addition, to explain the probable increase 
of frequency between the two orbits in the XMM-Newton observation of August 2002 (see 
§2.2), we have to assume that a small glitch of ~ 1 /(Hz occurred in the short 73 ks gap, 
which looks artificial. It is worth mentioning that the glitch hypothesis implies a smaller 
characteristic age (e.g., r ~ 60 kyr in the last example), bringing it closer to the SNR age, 
and larger magnetic field {B oc l/l^^^) and spin-down energy loss {E ex |/|). 
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3.2. Accretion as the source of the spin frequency variations? 

The spin evolution of a NS can also be affected by a flow of material streaming to the 
NS from an accretion disk. The analysis presented in §3.3 below rules out an accreting close 
binary, but the accreting material could be supplied from a "fossil" accretion disk formed 
from ejecta produced in the supernova explosion (see Marsden, Lingenfelter & Rotschild 
2001 for references). The accretion can proceed in two regimes (e.g., Prank, King, & Raine 
2002), depending on the relation between the corotation radius 

re = (G'M)i/3^27r/)-2/2 = 0.1 x 10^ cm, (4) 

(here / ~ 2.36 Hz is the spin frequency) and magnetospheric radius 

rm - 0.5(8 G M)-^/V^^^"^^^ - 0.9 x lO^rn'^^^ cm (5) 

(here m — lO^^rhu g s~^ is the accretion rate, ji — 10^°/i3o G cm^ is the NS magnetic 
moment, and the NS mass is assumed to be M = IAMq). If the accretion rate is so high 
that < Tc, the accreting matter can overcome the centrifugal barrier and reach the NS 
surface. In this "accretor regime" the torque exerted on the magnetosphere spins the NS up 
(/ > 0). At Tm > Tc, the centrifugal force at r = exceeds the gravitational force, so that 
accretion onto the NS surface is inhibited. In this "propeller regime" the infalling material 
is accelerated away from the magnetosphere reducing the angular momentum of the NS 
(/ < 0). At even lower m, approaches the light cyhnder radius, ric = c/ (27r/) = 2.0 x 10^ 
cm, where the pressure of the magneto-dipole radiation of the rotating NS takes over the 
accretion pressure and prevents any accretion ("ejector regime"). 

The torque caused by accretion in both the propeller and accretor regimes can be con- 
veniently approximated, at <^ ric, as 

K^ATTThriif^-f) (6) 

(e.g., Menou et al. 1999), where /eq = (27r)-^ {GMfl'^ r'^^''^ = 8.0 x lO-^m^f //^^^ Hz is the 
spin frequency at which the propeller- accretor transition occurs (/ = /eq at = r^. The 
corresponding frequency derivative, / = K{2'kT)~^ ^ is 

/ « 1.6 X 10-13 ^3/7^8/7 ^^^^ _ ^-2^ 

for the NS moment of inertia / = 10^^ g cm^. The curves of constant / values in different 
regimes are shown in Figure 4. 

If the observed variations of / are due to accretion, the accretion rate should vary with 
time. Moreover, to explain the very plausible spin-up between December 2001 and August 
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2002 (see §2.2), we have to assume, in the framework of this simple model, that 1E1207 
was in the accretor regime during at least part of that time interval. Figure 4 shows that in 
order to provide the corresponding / > +10~^^ s~^, the accretion rate should be high enough, 
rha > 10^^ g s~^ {TTia [< tji] IS the rate of accretion onto the NS surface). Because accretion 
onto the NS surface provides an additional luminosity AL ~ GMrhaR'^ — 1.9 x lO'^'^ma ie 
erg s~^, which substantially exceeds the observed X-ray luminosity, Lx — 1-5 x lO^^dl i 
erg s~^ in the 0.3-5 keV range {d2.i = (i/[2.1 kpc]), we have to assume that 1E1207 was 
in the propeller stage during all the available X-ray observations, while it perhaps was in 
the accretion stage (accompanied by an unnoticed outburst) between January and August 
of 2002. (This outburst plays the same role as the 7 //Hz glitch assumed in the glitching 
interpretation.) 

Assuming no other accretion outbursts occurred in the 3.45 yr of the Chandra and 
XMM-Newton observations, typical spin-down rates are in the range of |/| ~ 10~^^-10~^^ 
s""^. According to Figure 4, in the propeller regime they correspond to m ~ 10^^-10^^ 
g s~^ and n ~ 10^^-10'^° G cm^ (or B ~ 10^^-10^^ G for a centered magnetic dipole). 
The higher values of m from the above range are hardly plausible because such strong 
accretion onto the magnetosphere should be accompanied by emission (presumably, optically 
thin bremsstrahlung - see Wang & Robertson 1985) with a luminosity L ~ GMyhr^ ~ 
10^^ jjL^Q^'' erg s~^. The lack of such emission in the observed data means that m < 10^^- 
10^^ g s~^. At lower m (< 10^^ g s~^), the X-ray pulsar approaches the propeller-ejector 
boundary, where the propeller spin-down becomes less efficient but the radio pulsar can turn 
on. If m varies around this boundary, 1E1207 can occasionally manifest radio-pulsar activity. 

In the propeller regime, we can expect variations of the accretion rate comparable with 
its average value, to explain the variations of /. If a fraction of m accretes onto the NS 
surface (which is possible even in the propeller regime — e.g., Menou et al. 1999), it could 
contribute to the observed X-ray luminosity and lead to its variations. The fact that the 
observations of 1E1207 at different epochs do not reveal significant changes in its X-ray 
flux allows one to constrain this fraction. For instance, the difference between the fluxes 
derived from the two XMM-Newton observations does not exceed 1%, which corresponds to 
ALx < 1.5 X lO^^li erg s~^ and Ama ~ rha < 10^^ g s"^ -C m. Observations of 1E1207 
with different X-ray observatories show flux differences of up to 30%, but they can hkely be 
explained by systematic uncertainties in instrument calibrations. 

Strong constraints on the properties of the putative accretion disk can be obtained from 
optical observations of 1E1207. Using a standard approach (e.g., Frank et al. 2002), we can 
estimate the effective temperature of a geometrically thin, optically thick disk as 



T{r) = (3G'Mm)^/^(87rrV)-V^ = T,^{n.lrf'\ 



(8) 



- 11 - 



where 



Tin^l.SxlO^mlfVao K 



(9) 



is the temperature at the inner edge of disk, at Tin ~ r^. The spectrum of such a disk, in 
the blackbody approximation, is 



where i is the disk inchnation. We calculated the spectral flux and the corresponding V 
magnitude on a /j-ttl grid and plotted the lines of constant V in Figure 4 for cosi = 1, = 2.1 
kpc, adding a plausible extinction Ay = 0.7. The results are insensitive to the value of the 
outer radius rout when it exceeds ~ 10^-*^ cm (we used rout = lO"*^^ cm in this calculation). 

Figure 4 shows that, for the face-on disk, expected magnitudes in the accretor and 
propeller regimes arc V ~ 17-19 and V ~ 20-23, respectively. To reconcile the magnitude 
estimates with the reported limit, ^ > 25, we have to assume inclinations close to 90°, 
which looks rather artificial. On the other hand, the inclination should not be too close 
to 90° in order the disk periphery not to obscure the X-ray source: cosi > -ff(rout)/?'out ~ 
'i-0~'^ rhi^^^rl{^ i2i where H(r) oc r^/^ is the vertical scale- height of the disk, in the thin disk 



To conclude, the accretion hypothesis can, in principle, explain the observed variations 
of pulsation frequency and give some constraints on the magnetic field, but it implies that 
1E1207 is a transient X-ray source, which has not been observed, and requires extreme 
inclinations of the accretion disk. 



Another plausible explanation for the observed frequency variation is that 1E1207 resides 
in a binary system. The orbital motion would result in a periodic modulation of the observed 
frequency due to the Doppler shift. We can put some limits on the nature of the binary 
system and the companion of 1E1207 from the available X-ray and optical observations. 

Using our timing results, we can constrain the amplitude /„ of the orbital modulation 
and the orbital period Porb- The observed deviations from a steady spin-down imply 2 < 
fa ^ 20 //Hz (larger /„ values cannot be formally excluded but seem hardly plausible, given 
the range of the measured frequencies). To see a significant deviation from the steady spin- 
down over the 3.45 yr time span, the orbital period should not be larger than about 6 
years, for /„ < 20 /xHz. On the other hand, the orbital period should substantially exceed 




(10) 



model. 



3.3. A binary with a low-mass companion? 



- 12 - 



the time spans of the long observations of August 2002 (Tgpan = 3.4 d) and June 2003 
(^lipan = 9-3 d) that do not show frequency variations except for the hnear term, ft. To 
estimate a lower limit on Porb, we repeated the timing analysis of these two observations for 
a large number of models with sinusoidally modulated fit). The lower limit, naturally, grows 
with increasing /«, and it strongly depends on the assumed orbital phase. For example, we 
obtained Porb > (2 — 50) Tgpan for /„ = 2 /xHz. This gives the lowest possible period of 
^ 20 d, for a special choice of orbital phase, while the lower hmit is significantly larger, 
-Porb ^ 50-100 d, for a broad range of orbital phases. We also attempted to constrain the 
lower limit on Porb fitting the time dependence 

f{t)=fo + fot + fa Sin[27r (t - tp)/Porb] (11) 

(which assumes a circular orbit) to all five timing points shown in Figure 3 and requiring 
the model to be consistent with the timing data in each individual data set. This resulted 
in a 'global lower hmit' Porb ^ 60 d. 

Within the possible period boundaries, multiple combinations of orbital parameters 
can give "perfect fits" to the scarce timing observations. We present three representative 
solutions in Figure 5. These solutions correspond the following parameters in equation (11): 
/o = 65.4 /iHz, /o = -1.8 X IQ-^'^ s-^ = 4.8 /iHz, Po^b = 201.4 d and tp = -12.5 d (Model 
I), /o = 69.3 /iHz, /o = -8.0 x lO'^^ s^^, = 4.6 /xHz, P^rh = 937.8 d and tp = 18.8 d 
(Model II), and /„ = 65.6 /xHz, /o = -1.6 x 10"^^ s'^, /„ = -6.0 /iHz, Po^b = 595.2 d and 
tp = 111.0 d (Model HI), with time t counted from MJD 51500.0. Model I assumes that 
the correct frequency in the observation of June 2003 is given by peak A (see § 2.3 and Fig. 
1), whereas the more probable frequency of peak B is chosen for the two other models. We 
note that in these examples the characteristic pulsar age, Tc = — // (2/), ranges from 466 to 
2 335 kyr; it is still significantly larger than the 20-kyr upper limit on the age of the SNR. 
To reconcile the ages in the framework of this simple model, one would need a very large, 
and rather implausible, orbital modulation > 120 /iHz. 

The above estimates show that 1E1207 could be a wide binary with a long period, 
Porb ~ 0.2-6 yr, and the component separation, a ~ (0.3 — 3) (mi + 7712)^'^^ AU, where 
Ml = tjiiMq and M2 = m2MQ are masses of the NS and its companion, respectively. The 
plausible amplitudes of /„ correspond to rather small amplitudes of radial velocity, Vir = 
0.64 (/a/5 /iHz) km s^^. The mass function of the binary can be estimated as {1112 sini)^ (mi + 
m2)-2 = (P„rb/27rG') f i\ Mq 1 = 1.0 x 10-5(Porb/l yr)(/a/5 /iHz)^. For sini > 0.03, this 
gives 777,2 sinz = 0.027 (Porb/1 yr)^^^ (/a/5 /iHz) (7771/1. 4)^/^. We see that for reasonably large 
inclinations the secondary star should have a very low mass, in the range of brown dwarf 
or M-dwarf masses, with a lower limit as small as 6 Jupiter masses. The upper limit on 7712 
can be estimated from optical observations. For instance, the limiting magnitude V > 25 
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(Mereghetti et al. 1996) implies the absolute V magnitude My > 12.7 (for d — 2.1 kpc and 
a plausible extinction Ay = 0.7), corresponding to an M dwarf later than M3-M4 (mass 
777.2 ^ 0.25-0.3) or a white dwarf (of any allowed white dwarf mass, m2 ^ 1.4). 

Thus, we conclude that the timing results can be explained by the binary hypothesis 
if either the mass of the secondary is very low or the inclination is very small. Similar to 
the case of glitches, the intrinsic / (hence the characteristic age, spin-down energy loss, 
and magnetic field) can be quite different from the values inferred in Paper II from two 
Chandra observations. It should be noted that a low mass of the secondary implies that the 
pre-supcrnova binary had a very eccentric orbit, with the secondary close to aphelion at the 
moment of explosion, in order for the binary to survive. Therefore, such low-mass binaries in 
SNRs should be a very rare case. The only potential low-mass binary with a young NS is the 
central object of the SNR RCW 103, but hkely it is a short-period (Porb ~ 6.4 hr), accreting 
binary (Sanwal et al. 2002b). We note also that if 1E1207 is indeed in a wide binary system, 
we would not see manifestations of binarity other than the timing properties. 



4. Conclusions 

Our timing analysis of the Chandra and XMM-Newton observations, spread over a 3.45 
yr time span, has shown a non-monotonous frequency evolution of 1E1207. To explain the 
observed deviations from the steady spin-down, we have discussed three hypotheses. The first 
one, that the neutron star is a glitching pulsator, rcqiiircs an unusually large time-integrated 
amplitude of glitches to explain the observed variations in the spin frequency. The second 
hypothesis assumes variable accretion from a dim fossil (residual) disk around the pulsar and 
implies that 1E1207 is a transient X-ray source that can occasionally increase its brightness 
by at least three orders of magnitudes, which has never been observed. The deep limit on 
optical emission from 1E1207 imphes that such a disk is seen almost edge-on. The third 
hypothesis is that 1E1207 resides in a wide binary system with a low-mass companion that 
has not yet been detected. A binary with such a companion could survive the supernova 
explosion only at very special conditions. Thus, all three hypotheses imply rather exotic 
properties of 1E1207, but none of them can be firmly ruled out. We consider the binary 
hypothesis as somewhat more plausible than the other two, but only further observations can 
tell us which (if any) of these interpretations is correct. The most direct way to understand 
the true nature of 1E1207 is to monitor its timing behavior in a dedicated Chandra or XMM- 
Newton program. Additional useful constraints on the nature of the putative low-mass 
companion or accretion disk can be obtained from extremely deep optical/IR observations. 
Finally, it would be worthwhile to carry out a series of very deep radio observations to look 
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for a possible transient radio pulsar. 
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Fig. 1. — Probability density distributions p{f) for five observations of 1E1207 (/rcf = 
2.3577 Hz is the reference frequency). The insets in the two lower panels present the p{f) 
dependences in narrower frequency ranges. 
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Fig. 2. — Probability distributions p(A) for frequency difference A for two pairs of obser- 
vations: August 2002 - December 2001 (upper panel) and June 2003 - August 2002 (lower 
panel) . 
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Fig. 3. — Possible glitching scenario to explain the observed variations of the pulsation 
frequency (median frequencies with their 68% uncertainties are plotted) . The diamonds and 
triangles indicate the Chandra and XMM-Newton observations, respectively. The points A, 
B, and C are three timing solutions for the June 2003 observation (see Fig. 1). The sohd hne 
corresponds to / = —1.6 x 10^^^ s~^ between the glitches, and it assumes that the correct 
frequency in June 2003 is given by point B. The dots show the best straight-line fit with 
/ = —5.2 X 10~^^ s~^, assuming the correct frequency in June 2003 is given by point A. 
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Fig. 4. — Three regimes of interaction of a NS with accreting matter, for the NS spin 
frequency / = 2.36 Hz (see §3.2). The thin sohd curves are hues of constant / values 
(depicted near the curves, in units of s~^). They are obtained from equation (7) for the 
accretor and propeller regimes. In the ejector (radio-pulsar) regime, the frequency derivative 
is calculated as / = — (37r/)~^|U^rj^'^, assuming a magneto-dipole braking. The dotted parts of 
the curves in the propeller regime are drawn arbitrarily to demonstrate decreasing efficiency 
of propeller braking when the magnetosphere radius approaches the light-cylinder radius 
r\c. The long-dash curves correspond to constant V magnitudes of the predicted optical 
emission from an accretion disk (eq. [10]), for Ay — 0.7 and i — 0. 
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Fig. 5. — The same timing points as in Figure 3 are fitted with three models of binary 
motion (see §3.3 for the model parameters). The dashed lines show the steady spin-down 
components in the binary-model fits. The filled diamonds for the June 2003 observation 
indicate the frequency values chosen for each of the models at that epoch (see Fig. 1 and the 
text for details). 



